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Semimartingale reflecting Brownian motions (SRBMs) living in 
the closures of domains with piecewise smooth boundaries are of in- 
terest in applied probability because of their role as heavy traffic ap- 
proximations for some stochastic networks. In this paper, assuming 
certain conditions on the domains and directions of reflection, a per- 
turbation result, or invariance principle, for SRBMs is proved. This 
provides sufficient conditions for a process that satisfies the definition 
of an SRBM, except for small random perturbations in the defining 
conditions, to be close in distribution to an SRBM. A crucial ingredi- 
ent in the proof of this result is an oscillation inequality for solutions 
of a perturbed Skorokhod problem. We use the invariance principle 
to show weak existence of SRBMs under mild conditions. We also 
use the invariance principle, in conjunction with known uniqueness 
results for SRBMs, to give some sufficient conditions for validating 
approximations involving (i) SRBMs in convex polyhedrons with a 
constant reflection vector fleld on each face of the polyhedron, and 
(ii) SRBMs in bounded domains with piecewise smooth boundaries 
and possibly nonconstant reflection vector flelds on the boundary 
surfaces. 

1. Introduction. Semimartingale reflecting Brownian motions (SRBMs) 
living in the closures of domains with piecewise smooth boundaries are of 
interest in applied probability because of their role as heavy traffic diffusion 
approximations for some stochastic networks. The nonsmoothness of the 
boundary for such a domain, combined with discontinuities in the oblique 
directions of reflection at intersections of smooth boundary surfaces, present 
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challenges in the development of a rigorous theory of existence, uniqueness 
and approximation for such SRBMs. 

When the state space is an orthant and the direction of reflection is con- 
stant on each boundary face, a necessary and sufficient condition for weak 
existence and uniqueness of an SRBM is known [14]. This condition involves 
a so-called completely-5 condition on the matrix formed by the reflection 
directions. An invariance principle for such SRBMs was established in [15] 
and used in [16] to justify heavy traffic diffusion approximations for cer- 
tain open multiclass queueing networks. Loosely speaking, the invariance 
principle shows that, assuming uniqueness in law for the SRBM, a process 
satisfying the definition of the SRBM, except for small perturbations in the 
defining conditions, is close in distribution to the SRBM. 

For more general domains with piecewise smooth boundaries, some con- 
ditions for existence and uniqueness of SRBMs are known. In particular, for 
convex polyhedrons with a constant direction of reflection on each boundary 
face, necessary and sufficient conditions for weak existence and uniqueness of 
SRBMs are known for simple convex polyhedrons (where precisely d faces 
meet at each vertex in d-dimensions) and sufficient conditions are known 
for nonsimple convex polyhedrons, see [4]. For a bounded domain that can 
be represented as a finite intersection of domains, each of which has a C^- 
boundary and an associated uniformly Lipschitz continuous reflection vector 
field, sufficient conditions for strong existence and uniqueness were provided 
by Dupuis and Ishii [6]; in fact, these authors study stochastic differential 
equations with reflection which include SRBMs. Despite these existence and 
uniqueness results, a general invariance principle for SRBMs living in the 
closures of domains with piecewise smooth boundaries has not been proved 
to date. (We note that for the special case when the directions of reflection 
are normal, that is, perpendicular to the boundary, there are a number of 
perturbation results for reflecting Brownian motions. Our emphasis here is 
on treating a wide range of oblique reflection directions.) 

Motivated by its potential for use in approximating heavily loaded stochas- 
tic networks that are more general than open multiclass queueing networks, 
in this paper, we formulate and prove an invariance principle for SRBMs 
living in the closures of domains with piecewise smooth boundaries with 
possibly nonconstant directions of reflection on each of the smooth bound- 
ary surfaces. An application of the results of this paper to the analysis of 
an internet congestion control model can be found in [13] . An outline of the 
current paper is as follows. 

The deflnition of an SRBM and assumptions on the domains and direc- 
tions of reflection are given in Sections 2 and 3, respectively. Some sufficient 
conditions for these assumptions to hold are provided in Section 3. Section 
4 is devoted to proving the main result of this paper, namely, the invari- 
ance principle. A key element for our proof of this result is an oscillation 
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inequality for solutions of a perturbed Skorokhod problem; this inequality, 
which may be of independent interest, is proved in Section 4.1. In Section 5 
we give some applications of the invariance principle. We prove weak ex- 
istence of SRBMs under the conditions specified in Section 3. We also use 
the invariance principle, in conjunction with known uniqueness results for 
SRBMs, to give sufficient conditions for validating approximations involving 
(i) SRBMs in convex polyhedrons with a constant reflection vector field on 
each face of the polyhedron, and (ii) SRBMs in bounded domains with piece- 
wise smooth boundaries and possibly nonconstant reflection vector fields on 
the boundary surfaces. 

Beyond its possible use in justifying SRBM approximations for stochastic 
networks, the invariance principle might be used to justify numerical ap- 
proximations to SRBMs. A further possible extension of the results stated 
here would involve an invariance principle for stochastic differential equa- 
tions with reflection. The oscillation inequality for the perturbed Skorokhod 
problem and associated criteria for C-tightness described in Sections 4.1 
and 4.2 are likely to be useful for this. We have not developed such an ex- 
tension here as that would involve introduction of extra assumptions that 
would make the result less relevant for potential applications to stochastic 
networks. In particular, the approximating processes would involve stochas- 
tic integrals driven by a Brownian motion, whereas in stochastic network 
applications, the Brownian motion typically only appears in the limit. 

1.1. Notation, terminology and preliminaries. Let N denote the set of all 
positive integers, that is, N = {1, 2, . . .}, M denote the set of real numbers, 
which is also denoted by (—00,00), denote the nonnegative half-line, 
which is also denoted by [0,oo). For x € M, we write for the absolute 
value of X, [x] for the largest integer less than or equal to x, for the 
positive part of x. For any positive integer d, we let M*^ denote d-dimensional 
Euclidean space, where any element in is denoted by a column vector. 
Let II • II denote the Euclidean norm on W^, that is, ||rc|| = (X)f=i )^''^ 
X G M.'^, and (•, •) denote the inner product on R^, that is, {x, y) = X]f=i Xiyi, 
for x,y € M'^. We note that for any x € W^, \\x\\ < J2i=i \xi\. Let R^^, denote 
the positive orthant in M'^, that is, R'l = {x e R'^ : Xi > 0,1 < i < d}. Let 
13{S) denote the Borel cj-algebra on 5 C M'^, that is, the collection formed 
by intersecting all Borel sets in M'^ with S. Let dist(x, 5) denote the distance 
between x G M'^ and 5 C M"*, that is, dist(x, S) = inf{|| x — y\\:yG S}, with the 
convention that dist(x,0) = 00 for x G M'^. Let Ur{S) denote the closed set 
{x G M"^ : dist(x, S) < r} for any r > and S C R'^, where if 5 = 0, Ur{S) = 
for all r > 0. Let Br{x) denote the closed ball {y G R'^ : ||y — x|| < r} for any 
X G M'^ and r > 0. For any set S C M*^, we write 5 for the closure of S, S° for 
the interior of S and dS = S\S° . For a finite set S, \S\ denotes the number 
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of elements in S. For any v G W^, v' denotes the transpose of v. Inequalities 
between vectors in should be interpreted componentwise, that is, if u, G 
W^, then u < {<)v means that Ui < {<)vi for each i G {!,..., d}. For any 
matrix A, let A' denote the transpose of A. For any function x:M+ W^, 
x{t—) denotes the left limit of x at t > whenever x has a left limit at t; 
unless explicitly stated otherwise, x(0— ) = 0, where is the zero vector in M'^. 
For any function x : R+ M*^, we let Ax{t) = x(t) — x(t—) for t € M4- when 
x(t—) exists. We let be the constant deterministic function 3;:M+ ^ M'^ 
such that x{t) = for all t € M+. 

A domain in M"^ is an open connected subset of W^. For each continuously 
differentiable function / defined on some nonempty domain S" C M'^, V/(x) 
is the gradient of / at x G S". For each x G M*^, a neighborhood 14 of x is a 
bounded domain in that contains x. For any nonempty domain 5 C M*^, 
we say that the boundary dS of S is C^, if for each x G dS there exists a 
Euclidean coordinate system Cx for M*^ centered at x, an r^; > 0, and a once 
continuously differentiable function ipx iM'^"^ —>■ M such that (px{0) = and 

SnBr^{x) = {z = {Zi,.. .^ZdY in Cx'-Zd > ^x{zi,.. . , Zd^i)} (1 Br^{x). 

Then, for x G dS, the inward unit normal to dS at z G dS PI B,.^{x) is given 
in the coordinate system Cx by 

n{z) = — - — {-Vipx{zi,...,Zd-iy,iy, 

(1 + \\\/ipx{Zl,...,Zd^lW)'-'^ 

where \/ipx{zi, . . .,Zd-i) = (^, • • • , ^^y{zi, . . .,Zd-i). For any nonempty 

convex set 5 C ffi'^, we call a vector n G M'^\ {0} an inward unit normal vector 
to 5 at X G dS if ||n|| = 1 and {n,y — x) >0 for all y G 5. Note that such a 
vector need not be unique. 

All stochastic processes used in this paper will be assumed to have paths 
that are right continuous with finite left limits (abbreviated henceforth as 
r.c.1.1.). A process is called continuous if almost surely its sample paths 
are continuous. We denote by D([0, 00), M*^) the space of r.c.1.1. functions 
from [0,00) into R'^ and we endow this space with the usual Skorokhod 
Ji-topology (cf. Chapter 3 of [7]). We denote by C([0, oo),M'^) the space 
of continuous functions from [0, 00) into W^. The Borel a-algebra on either 
L'([0,oo),M'^) or C([0,cx)),M'^) will be denoted by M'^. The abbreviation 
u.o.c. will stand for uniformly on compacts and will be used to indicate 
that a sequence of functions in D([0, 00), M"^) (or C([0, oo),M'^)) is converg- 
ing uniformly on compact time intervals to a limit in Z)([0, oo),R'^) (or 
C([0, oo),]R'^)). Consider W^,W^,...,W, each of which is a d-dimensional 
process (possibly defined on different probability spaces). The sequence 
is said to be tight if the probability measures induced by the 
on the measurable space {D{[0,oo),R'^),M'^) form a tight sequence, 
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that is, they form a weakly relatively compact sequence in the space of 
probability measures on {D{[0,oo),M'^),M'^). The notation 'W" W" wih 
mean that, as n — > cxd, the sequence of probability measures induced on 
(-D([0, oo), M'^), Al*^) by {M^""} converges weakly to the probability measure 
induced on the same space by W. We shall describe this in words by saying 
that converges weakly (or in distribution) to W as n ^ oo. The se- 
quence of processes {VF"}^]^ is called C-tight if it is tight, and if each weak 
limit point, obtained as a weak limit along a subsequence, almost surely has 
sample paths in C([0, cx)),M"'). The following proposition provides a useful 
criterion for checking C-tightness. 

Proposition 1.1. Suppose that, for each n£N, is a d- dimensional 
process defined on the probability space (r2",jF",i-'"). The sequence {W^}'^^i 
is C-tight if and only if the following two conditions hold: 

(i) For each rj > and T >0, there exists a finite constant M^i^t > 
such that 



(1) liminfP"! sup <M„t1 > 1 

lo<t<T J 

(ii) For each e > 0, t] > and T > 0, there exists A G (0,T) such that 

(2) limsupP"{w;r(T^",A) >e} <r/, 

n— >oo 

where for x e L'([0, cx3),M'^), 

(3) wt{x, X) = supl sup \\x{u) — x{v)\\ -.0 <t < t + X <T > . 

iu,vG[t,t+X] J 

Proof. See Proposition VI.3. 26 in [12]. □ 

A d-dimensional process W is said to be locally of bounded variation if 
all sample paths of W are of bounded variation on each finite time interval. 
For such a process W, we define V(VF) = {V{W){t),t > 0} such that for each 
t >0, 

v{w){t) = \\wm 

+ supj^ \\Witi) - W{t^-i)\\ ■.0 = to<ti<---<ti = t,l>l^. 

A triple {Q,T,{J^t,t > 0}) will be called a filtered space if $7 is a set, 
^ is a cr-algebra of subsets of $7, and {J-t,t > 0} is an increasing family of 
sub-a-algebras of that is, a filtration. Henceforth, the filtration {J-t,t > 0} 
will be simply written as {Tt}. If P is a probability measure on {Vt,T), then 
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(fi, J- ^ {-^t}, P) is called a filtered probability space. A d-dimensional process 
X = {X{t),t > 0} defined on {^l,J^,P) is called {.7^t}-adapted if for each 
t > 0, X{t) -.Q is measurable when Q is endowed with the cr-algebra 

Given a filtered probability space {Tt},P), a vector fieR'^,adxd 

symmetric, strictly positive definite matrix T, and a probability distribution 
u on (M'^, an {^t}-Brownian motion with drift vector /i, covariance 

matrix T, and initial distribution u, is a d-dimensional {^Ftj-adapted process 
defined on {0,,J^,{J^t},P) such that the following hold under P: 

(a) X is a d-dimensional Brownian motion whose sample paths are almost 
surely continuous and that has initial distribution u, 

(b) {Xi{t) — Xi{0) — fiit, J^ti i > 0} is a martingale for i = 1, . . . , d, and 

(c) {{Xi{t) - Xi{<d) - tiit){Xj{t) - Xj{Q) - iijt) - Tijt,Tut > 0} is a mar- 
tingale for i, j = 1, . . . , d. 

In this definition, the filtration {Tt} may be larger than the one generated 
by X\ however, for each f > 0, under P, the cr-algebra J-t is independent of 
the increments of X from t onward. The latter follows from the martingale 
properties of X. li u = 5x, the unit mass at x € M'^, we say that X starts 
from X. 

2. Definition of an SRBM. Let G = PIigi Gi be a nonempty domain in 
W^, where T is a nonempty finite index set and for each i ^ I, Gi is a 
nonempty domain in W^. For simplicity, we assume that I = {1,2,...,!} 
and then \I\ = I. For each i G T, let 7*(-) be a vector valued function defined 
from into W^. Fix S M*^, T a d x d symmetric and strictly positive 
definite covariance matrix and i' a probability measure on (G, B{G)), where 
B{G) denotes the cr-algebra of Borel subsets of the closure G of G. 

Definition 2.1 (Semimartingale reflecting Brownian motion). A semi- 
martingale refiecting Brownian motion (abbreviated as SRBM) associated 
with the data (G, /x, F, {7*, i gT},u) is an {J^t}-adapted, d-dimensional pro- 
cess W defined on some filtered probability space {i^,J^, {J-'t},P) such that: 

(i) P-a.s., Wit) = X{t) + /(o,i] lKW{s)) dYi{s) for all t > 0, 

(ii) P-a.s., W has continuous paths and W{t) G G for all t > 0, 

(iii) under P, X is a d-dimensional {jF^j-Brownian motion with drift 
vector /i, covariance matrix F and initial distribution v, 

(iv) for each i £T, Yi is an {.7^t}-adapted, one-dimensional process such 
that P-a.s., 

(a) YiiO) = 0, 

(b) Yi is continuous and nondecreasing. 
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(c) Yi{t) = /(o J] l{wis)edG,ndG} dMs) for all t > 0. 

We shall often refer toY = {Yi,i £ 1} as the "pushing process" associated 
with the SRBM W. When v = 6x, we may alternatively say that W is an 
SRBM associated with the data (G, /i, F, {7*, i € I}) that starts from x. We 
will call {W,X,Y) satisfying Definition 2.1 an extended SRBM associated 
with the data (G, n, T, {7*, i € I}, u). 

Loosely speaking, an SRBM behaves like a Brownian motion in the inte- 
rior of the domain G and it is confined to G by instantaneous "reflection" (or 
"pushing") at the boundary, where the allowed directions of "reflection" at 
X G dG are convex combinations of the vectors 7* (2;) for i such that x G dGi. 
Under the assumptions imposed on G and {7*,? S T} in Sections 3.1 and 3.2 
below, at each point on the boundary of G there is an allowed direction of 
reflection that can be used there which "points into the interior of G." We 
end this section by introducing a related set- valued function Z{-) and show 
a key property of it. 

Definition 2.2. For each 2; G R'^, let l{x) = {iel:x£ dGi}. 

The set-valued function X(-) has the following property called upper semi- 
continuity on dG. 

Lemma 2.1. For each x G dG, there is an open neighborhood Vx of x in 
R'^ such that 

(4) I{y)cl{x) forallyeVx. 

Proof. We prove this lemma by contradiction. Suppose that the func- 
tion Z(-) does not satisfy (4). Then there is a point x G dG such that there 
is no open neighborhood Vx oi x such that 2^(y) C I{x) for all y G T4. Since 
the index set I is finite, there is an index /c G X \ I{x) and a sequence of 
points {Un} C M"^ such that ||y„, — < ^ and k G for each n > 1. Hence 

Un G dGk for all n > 1. Since dGk is closed and yn ^ x as 00, we con- 
clude that X G dGk. This implies that k G I{x), which is a contradiction, as 
desired. □ 

3. Assumptions on the domain G and the reflection vector fields {7^^}. 

3.1. Assumptions on the domain G. We henceforth assume that the 
domain G satisfies assumptions (Al)-(A3) below. In the case when G is 
bounded, assumptions (A2)-(A3) follow from assumption (Al) (see Lem- 
mas A.l and A. 2 in the Appendix for details). If the domain G is a convex 
polyhedron satisfying assumption (Al), then assumptions (A2)-(A3) hold 
by Lemma A. 3 in the Appendix. 
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(Al) G is a nonempty domain in R"^ with representation 

(5) G=f]G„ 

iei 

where for each i £2, Gi is a nonempty domain, Gj ^ Mf^, and the boundary 
dOi of Gi is C^. For each i G I, we let be the unit normal vector field 
on dGi that points into Gi. 

(A2) For each e S (0, 1) there exists R{£) > such that for each i £l, 
X G dGi n dG and y € G satisfying [[x — ?/|| < R{e), we have 

(6) (n^x),!/ — x) > —e\\x — y\\. 

(A3) The function D : [0, oo) [0, c5o] defined such that D{0) = and 

(7) D{r) = sup sup J dist ( x, C] {dGj ndG)]:x£ fl Ur{dGj n dG) \ 
for r > 0, satisfies 

(8) D{r)^0 asr^O. 



Remark. Assumption (A2) is reminiscent of the uniform exterior cone 
condition (cf. [9], page 195). We say that a region G CW^ satisfies a uni- 
form exterior cone condition if for each xq G dG, there is a truncated closed 
right circular cone Vxq, with nonempty interior and vertex xq, satisfying 
Vxg nG = {xq}, and the truncated cones Vxq are all congruent to some fixed 
truncated closed right circular cone V. By comparing assumption (A2) with 
the uniform exterior cone condition, we see that assumption (A2) implies 
the uniform exterior cone condition. On the other hand, under assumption 
(Al), assumption (A2) is implied by a family of uniform exterior cone condi- 
tions where for each e £ (0, 1), the axis of the truncated closed right circular 
cone at X G dG is along the vector — n*(x) and all of the truncated closed 
right circular cones are congruent to a truncated closed right circular cone 
whose height and base radius are R{£) and R{e){-p — 1)^/^ respectively. 
Assumption (A2) holds automatically if G is convex. We also note that as- 
sumption (A2) is strictly weaker than the uniform exterior sphere condition. 
The definition of the uniform exterior sphere condition is similar to that for 
the uniform exterior cone condition where a closed ball with xq on its bound- 
ary takes the place of the truncated closed right circular cone Vx^- It can 
be checked that for the domain G = {{x,y) G :?/ < with a G (1,2), 

the uniform exterior sphere condition fails to hold, but assumption (A2) 
holds. In fact, at the point (0, 0) G M^, there is no r > and y G such that 
Br{y)ndG = {{0,0)}. 
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Remark. For the definition of D{-) in (A3), we adopt the convention 
that the supremum over an empty set is zero and dist(a;, 0) = oo. Since dGi PI 
dG 7^ for at least one i(zl, the function D{-) satisfies linir^oo D{r) = oo. 
Furthermore, D{ri) < D{r2) whenever ri,r2 G [0, oo) and ri < r2. Assump- 
tion (A3) requires that for any nonempty subset J <Zl, the intersection of 
tubular neighborhoods of the boundaries {dGj H dG : j S J^} given by the 
set C\j^jUr{dGj n dG) "converges" to the intersection of the boundaries 
given by the set CljejidGj n dG) as r approaches 0. Property (8) need 

not always hold. For example, let Gi = {{x,y) £M? :y < e~^^/^, x G M} and 
G2 = {ix,y) eR'^:y>0,x£R}. Then dGi n dG2 = 0. But for each r > 0, 
Ur{dGi) n Ur{dG2) / 0. Hence D{r) = 00 for each r > 0. 

3.2. Assumptions on the reflection vector fields {7*}. We henceforth as- 
sume that there are vector fields {Yi')^''- satisfying assumptions (A4)- 
(A5) below. 

(A4) There is a constant L > such that for each i ^I, 7*(-) is a uni- 
formly Lipschitz continuous function from into with Lipschitz con- 
stant L and ||7'(x)|| = 1 for each x € W^. 

(A5) There is a constant a £ (0,1), and vector valued functions 6(-) = 
(6i(-), . . . , fei(-)) and c(-) = (ci(-), . . . , ci(-)) from dG into such that for 
each X G dG, 

(i) EiGX(x)^i(a;) = 1> 

(9) min ( hi{x)n\x),^^ {x)) > a, 

(ii) EiGJ(x)Ci(x) = 1, 

(10) min ( Ci(x)7*(2;),re-^(rc) ) > a. 

We note here for future use that by (A4), if we set Po = then for 
any x,y € satisfying ||x — y\\ < po, we have ||7*(x) — 7*(y)|| < a/4 for 
each i G Z. So for each < p < pq/A, by (9)-(10) and the normalization of 
6(-),c(-),7*(-),n-'(-) for i,jGT, we obtain 

(11) Hr< ■'^^J'^ 1 Hx)n\x),^i{y))>a/2 



\iel{x) 



and 



(12) ijif, ^( E Ci{x)-f\y),n^{x))>a/2. 

x€dGjeI{x)yeBip{x)\.^^^ I 
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The use of Bip{x) here is related to the form in which this is used in Section 
4.1. 

Remark. Assumption (A4) is equivalent to (3.4) in [6] when G is bounded. 
Property (10) means that, at each point x G 9G, there is a convex combi- 
nation 7(x) = Z]iGX(a;) Ci(ic)7*(x) of the vectors {'y^{x),i ^T{x)} that can 
be used there such that 7(x) "points into" G. Property (9) is in a sense a 
dual condition to property (10), where the roles of 7' and n* are reversed 
for i £l{x). This property (9) is used in showing the oscillation inequality 
in Theorem 4.1 below. Assumption (A5) is an analogue of Assumption 1.1 
in [4]. When G is bounded, (10) is similar to condition (3.6) in [6] (we as- 
sume some additional uniformity through the lack of dependence of a on 
x). 

It is straightforward to see using the triangle inequality that the following 
condition (A5)' implies (A5). 

(A5)' There is a € (0, 1) and vector valued functions 5, c from dG into 
such that for each x E dG, 

(i) J2iej{x) bi{x) = 1, and for each i G T{x), 

(13) b,ix){n\x),y{x))>a+ ^ bjix)\{n^ ix),Yix))\, 

j€lix)\{i} 

(ii) J2iej{x) (^i{x) = 1, and for each i G T{x), 

(14) Ci{x){Y{x),n\x))>a+ Cj{x)\{-f^ {x),n\x))\. 

jei{x)\{i} 

Condition (A5)'(ii) is similar to condition (3.8) in [6], although here we 
assume additional uniformity through the lack of dependence of a on 2; . As 
noted in [6], their condition (3.8) can be phrased in terms of a nonsingu- 
lar M-matrix requirement [2]. (This is sometimes also called a generalized 
Harrison-Reiman type of condition [10].) Since that nonsingular M-matrix 
property is invariant under transpose, and this property for the transpose 
corresponds to a local form of (A5)'(i), one might conjecture that there is an 
equivalence between the existence of a nonnegative vector valued function b 
such that (A5)'(i) holds for each x € dG and the existence of a nonnegative 
vector valued function c such that (A5)'(ii) holds for each x S dG. Indeed 
we have the following lemma. We have stated the two (equivalent) condi- 
tions (i) and (ii) in specifying (A5)' to preserve a parallel with (A5) and 
since both properties can be useful in proofs. Furthermore, in light of the 
following lemma, verifying either condition suffices for both to hold. 
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Lemma 3.1. There is a constant a € (0, 1) and a vector valued function 
b:dG M?j_ such that (A5)'(i) holds for each x G dG if and only if there 
is a constant a € (0, 1) and a vector valued function c: dG —>■ such that 
(A5)'(ii) holds for each x € dG. 

Proof. We just prove the "if" part; the "only if" part can be proved 
in a similar manner. 

We suppose that there is a constant a G (0, 1) and a vector valued function 
c:dG such that (A5)'(ii) holds for each x € dG. For fixed x € dG, 

consider the square matrix A{x) whose diagonal entries are given by the 
(positive) elements (n*(x), 7*(x)) for i sT(x) and whose off-diagonal entries 
are given by —\{n'^{x),^^{x))\ for i G Z{x),j € 1{x), j ^ i. Let E be the square 
matrix having the same dimensions as A(x) and whose entries are all equal 
to one. By the theory of M-matrices (see [2], Chapter 6, especially condition 
(M35)), condition (ii) of (A5)' implies that A{x) — is a nonsingular M- 
matrix, that is, A{x) — ^E has nonnegative diagonal entries and nonpositive 
off-diagonal entries and it can be written in the form s{x)I — B{x) where 
B{x) is a matrix with nonnegative entries and s{x) > is a constant that is 
strictly larger than the spectral radius of B{x). 

Since the nonsingular M-matrix property is invariant under transpose (cf. 
(G21) in Chapter 6 of [2]), then A'{x) — ^E is also a nonsingular M-matrix. 
Hence, there is a vector b{x) = {bi{x) : i G T{x)) with nonnegative entries such 
that {A'{x) - ^E)b{x) > (cf. (I27) in Chapter 6 of [2]). We can extend b{x) 
to an I-dimensional vector 5(x) and normalize it so that J2i£2[x) hi^) = 1- 
Then (A5)'(i) holds with | in place of a. □ 

4. Invariance principle. In this section we state and prove an invariance 
principle for an SRBM living in the closure of a domain G with piecewise 
smooth boundary and having associated reflection fields {7*,* S^T}, where 
G, {7*,zGT} satisfy assumptions (Al)-(A5) of Section 3. (These assump- 
tions hold throughout this section.) We shall first state a preliminary result 
called an oscillation inequality (see Theorem 4.1), then we use it to prove 
a tightness result (see Theorem 4.2). Finally, we establish the invariance 
principle (see Theorem 4.3). 

4.1. Oscillation inequality. The following oscillation inequality is the key 
to the proof of the tightness result claimed in Theorem 4.2. In this subsec- 
tion, for any < ti < t2 < 00 and any integer k>l, D{[ti,t2],M.^) denotes 
the set of functions w : [ti,t2] that are right continuous on [ti,t2) and 

have finite left limits on (ti,t2]- For w G D{[ti,t2],^^), 

(15) Osc(u', [ti,t2]) = sup{||u'(t) - w{s)\\:ti <s<t< ^2}, 

(16) Osc{w, [ti,t2)) = sup{\\w{t) - w{s)\\:ti <s<t< t2}. 
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Note that we do not explicitly indicate the dependence on k in the notation. 

Recall the constants a,L from assumptions (A4)-(A5), the functions R{-) 
from assumption (A2) and D{-) from (7). Let po = 

Theorem 4.1 (Oscillation inequality). There exists a nondecreasing func- 
tion n : (0, oo) (0, oo] satisfying Il{u) ^0 as u^O, such that H depends 
only on the constants I, a and the function D{-), and such that whenever 

< p < min{f ,^^^}, 0<5<f, 0<s<t<cx), w,x ^ D{[s,t\,W^) and 
y G D([s,t],]R^) satisfy: 

(i) w{u) G Bp{xo) n Us{G) for all u G [s,t], for some xq G G, 

(ii) w{u) = w{s) + x{u) - x{s) + Eiex/{s,«]7*(^(^'))%(^^) for all u G 

(iii) for each i G X, 

(a) yi{s)>0, 

(b) yi is nondecreasing and Ayi{u) < 5 for all u G {s,t], 

(c) =yi(s) l{^„(^,)eC/^{9G,nc»G)}'^?/^(^^) for all u G 

(iv) L'(n(0sc(x,[s,t]) + (5))<f, 
i/ien toe /lave that the following hold: 

(17) Osc(w;, [s,t]) <n(Osc(x, [s,t]) + (5), 

(18) Osc(?/, [s, t]) < n(Osc(x, [s, t]) + (5). 

Proof. Let 

no(n) = n for all u> 0. 
Define Ilm ■ (0, oo) ^ (0, oo], m = 1, . . . ,1, inductively such that 

Umiu) = n^-iiu) + (I + 2)u + (l + ^) {D{n^_i{u) + (I + 2)n) + 2^). 

Here the sum of any element of [0, oo) with oo is oo and D{oo) is defined 
to equal oo. For each m = 0, 1, ... ,1, the function Ilm is nondecreasing and 
depends only on I, a and D{-). For each m = 1, . . . ,1 and u> 0, 11^-1 (it) < 
Umiu). By assumption (A3), we conclude (using an induction proof) that 

^m{u) -^0 as n ^ 0, for m = 0, 1, . . . , I. 

Let n(-) = ni(-). 

Fix < p < min{^, ^^^}, 0<5<f, 0<s<t<oo. Suppose that 
w,x G L>([s,t],M"') and y G L>([s,t],]R^) satisfy (i)-(iv) in the statement of 
Theorem 4.1. For each nonempty interval [ti,t2] C [s,t], let 

■^[tiM = {i -wiu) ^ Us{dGi n dG) for some u G [ti, i2]}) 
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the indices of the boundary surfaces that w{-) comes close to in the time 
interval [ii,i2]- For each < m < I, define %n = {[ti,t2] C [s,t] : ^^jl ^ 
m}. Note that under the partial ordering of set inclusion, %n. increases with 
m. To prove the theorem, we will prove by induction that for each < m, < I 
and each interval [ti,t2] G %n, (17)-(18) hold with [ti,t2] in place of [s,t] and 
Ilm{-) in place of n(-). The result for m = I yields the theorem. 

Suppose that m = 0. Then Tq = {[ii,t2] C [s,t] : l^ft^^tj] | = 0}. Fix an inter- 
val [ii,t2] S %■ Since I[ti,t2] — ^ ^nd (iii)(c) holds, the function y does not 
increase on the time interval {ti,t2], that is, yi{t2) — yi{ti) = for all i^I. 
Then, for ti < u < f < t2, 

(19) w{v) — w{u) = x{y) — x{u). 
So in this case, 

(20) Osc('W,[ti,t2]) =0sc(x,[ti,t2]) <Osc(a::,[ti,t2]) + 5, 

(21) Osc(?/, [tiM]) = < Osc(rE, [tiM) + ^■ 

Thus, (17)-(18) hold with no(-) in place of n(-) and [ti,t2] in place of [s,t] 
for each interval [ti,t2] G Tq. 

For the induction step, let 1 < m < I and suppose that (17)-(18) hold 
with Hm-i{') in place of n(-) and [ti,t2] in place of [s,t] for each interval 

Now fix [ti,t2] £ Tm- If |2^[ti,t2]l < "I — 1, then [ti,i2] G Tm-i and so by 
the induction assumption we have that (17)-(18) hold with [ii,t2] in place 
of [s,t] and Iirn~i{-) [and hence !!„,(•)] in place of n(-). Thus, it suffices 
to consider [ti,t2] C [s,t] such that llj^^ ^j]! = "i. For i ^1[ti,t2]^ by (iii)(c), 
yi{t2) — yi{ti) = 0, and so by (ii), for ti < u < ?; < t2i we have 




Let Yl.m{u) = Iim-i{u) + (I + 2)ti for all u > 0, and rj = Osc(x, [ti,t2]) + 5- 
For any M € (0, oo] and any nonempty set J <ZZ, let 

= G M-^* : dist(z, dGi n dG) < M for all i € J}. 

Note that Fj = when there is an i € ^ such that dGi H dG = 0. Since 
nm(') < !!„,(•) < n(-), D{-) and n(-) are nondecr easing, and Osc(x, [ti,t2]) < 
Osc(x, [s,t]), we have by (iv) that 

(23) D(jiM) < DiUM) < DiUiij)) < ^. 

Note that this implies Umiv) < oo since D{oo) = oo. 
We now consider two cases. 
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Case 1. Suppose that w{r) G Fj^^^^^^ for all r E [ti,t2] 
Fix u, V such that ti < u < v <t2- Since we have that 



by the definition of D{-) and (23), there is z G Cljex^^ ^ ^i'^^j ^ ^^^h 
that 

(24) \\w{v) - z\\ < DiUM) < ^. 

For each r £ [ti,t2], by (i) we have that 'w{r) € Us{G), and so there is z** E G 
such that 

\]w{r) - z''\\<2S. 
Hence by (i) and (24) we have 

\\z'^ — z\\ < Wz"^ — tt;(r)|| + \\w{r) — xo\\ + ||xo — u^('v)|| + \\w{v) — z\\ 

(25) 

<26 + p + p + p/2<4p< R{a/4:) 

and 

\\w{r) — z\\ < \\w{r) — xq\\ + ||a;o — tf(i')|| + Wwiv) — z\\ 

(26) 

<P + P + P/2<4:P. 

By (6) and (25) we have 

(27) {n^ (z), z — z^) < —\\z — z^W for each j E X(z) and r E [ti, t2]- 

Note that I{z) dXjj^^^^]- Recalling the definition of b(-) from assumption 
(A5), on dotting the vector ^jizx{z)^'j{z)n^ [z) with both sides of (22) and 
rearranging, we obtain 

I (T. b,{z)n^{z),f{w{r)))dy,{r) 

(28) = b,iz){n^{z),w{v)-w{u)) 

- Y bj{z){n^z),x{v) -x{u)). 

So by (11), (22), (24)-(28), and the fact that Ejex(^) bj{z) = 1, hj{z) > for 
j E I, we have 
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< E / ( E h,{z)n\z),^\w{r)))dyi{r) 

= h^{z){n^{z)Mv)-z)+ h,{z){n\z),z-z^) 

+ h{z){n^{z),z^-w{u))- Y hj{z){n\z),x{v)-x{u)) 
jeiiz) jeiiz) 

< D{UM) + j\\z- +26 + \\x{v) - x{u)\\ 
<D{I[M)+26+\\x{v)-x{u)\\ 

+ ^{\\z - w{v)\\ + \\w{v) - w{u)\\ + \\w{u) - z"| 
<D{UM)+26+\\x{v)-x{u 



+ ^(d{UM) + Mv)-x{u)\\+ Y {y^{v)-yi{n)) + 26\ 

<fi + ^){D(iiM) + '^s + Mv)-x{u)\\} + j Y imiv)-m{u)). 



Hence 



^ Y {y^{v)-y,{u))<{l + ^){D{n^{l^)) + 25+\\x{v)-x{u)\\} 



4 -cT 



<(l + ^){D(n„(7?))+2r/}. 



4 

On multiplying through by ^, we obtain 

(29) Y (yK^^) - y^{u)) ^ (i + ^) + 2r?} < n„(7?). 

Hence, by (29) and the fact that for any x G W^, \\x\\ < J2i=i \xi\, we have 

(30) Osc{y,[ti,t2]) <nm{0sc{x,[ti,t2]) + 5), 

and by (22), (29) and the definitions of Hm(-) and Ilm{-), we have 

Osc(u;, [h,t2]) < Osc(2;, [h,t2]) + f 1 + {£'(n^(r/)) + 2r?} 



a J 

<n„(Osc(x,[ti,t2]) + 5), 

as desired. 
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Case 2. Suppose that there is ts G [ti,t2] such that w{t-i) ^ i^^'"^''^ 

Define a = inf{n G [ti,t2\:w{u) ^ Fj^^^f^}- Then a < ^2- For each u G 

[ti,(T), w^u) G Fj^^^^ and so by a similar analysis to that for Case 1, we 
obtain for each v £[ti,a), 

Osc{w, [h,v]) <r]+(^l + ^^ {D{Tim{v)) + 2ri) 



and 



Osc(y, [h,v]) < (1 + ^) iD{UM) + 2??). 

By the right continuity of paths we have w{a) ^ Fj^'^^^K Then there is an 

i GT[j^ j2] such that dist (it; (cr), dGi D dG) > Ilm{ri), and it follows that w 
does not reach Us{dGi (1 dG) during the interval [a, t2]- To see this, let 
r = inf{n G [o",t2] : dist(w;(u), 9Gi dG) < 5} with the convention that the 
infimum of an empty set is 00. If t < t2, then by the right continuity of 
w{-) and since Ylmiv) > we have t > a and dist(u;(T), SGj fl dG) < 6. Also, 
since I^T^j^^^^] | = ''ti-, we have [a, u] G %n-i for each u G [a, r). By the induction 
assumption and letting u ^ r, we have ||i(;(t— ) — < JIm-i{ri). By (ii), 

(iii)(b) and since ||7*(-)|| = 1, we have 

||A«;(r)|| < i|Ax(T)|| + ^ Ay,(r) < Osc(x, [hM) + IS< Irj. 
Then simple manipulations yield 

dist{w{a),dGir\dG) < ||u;(cj) - u;(r-)|| + \\Aw{t)\\ +dist{w{T),^G^n^G) 
<n„„i(?7) +Ir/ + 5 

<n^(r/). 

This contradicts the fact that dist{w{a), dGiCidG) > ILm{r]), and so confirms 
that w does not reach Us{dGi HdG) in [(7,^2] • Thus we must have [ex, ^2] £ 
%n~i. Hence we have by the induction assumption that 

Osc{w,[ti,t2]) < sup Osc{w,[ti,v]) + \\Aw{a)\\ + Osc{w,[a,t2]) 
<r]+(l + -] (DiUM) + 2??) + Ir/ + H^-i (r/) 



a, 

<n„(0sc(x,[ti,t2]) + <^) 
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and 

Osc{y,[ti,t2]) < sup 0sc{y,[ti,v]) + \\Ay{a)\\+0sc{y,[a,t2]) 

< (^1 + 1^ iD{Timiv)) + 2?7) + Ir? + n^-i(??) 

<n„,(Osc(x,[ti,t2]) + <5). 
On combining all of the cases above, we have 

(31) Osc(u;, [h,t2]) < n^(Osc(x, [h,t2])+6), 

(32) Osc(y, [tiM) < n^(Osc(x, [ti,t2])+6). 
This completes the induction step. □ 

Remark. The proof of the above theorem was inspired by the proof of 
Lemma 4.3 of [4]. Because of the condition (i) in Theorem 4.1, the oscilla- 
tion inequality given here is localized. Similar, but nonlocalized, oscillation 
inequalities were proved in [15] when G = and in [3] for a sequence of 
convex polyhedrons; in these cases, the direction of reflection was constant 
on each boundary face. 

4.2. C -tightness result. Throughout this subsection and the next, we 
suppose that the following assumption holds in addition to (Al)-(A5). 

Assumption 4.1. There is a sequence of strictly positive constants {6^}'^^i 
such that for each positive integer n, there are processes W^^W^^X^^a^ 
having paths in Z)([0, oo),M"') and processes Y^,Y'^,j3'^ having paths in 
Z)([0, cxd),]R^) defined on some probability space {Q'^,J^,P"') such that: 

(i) P"-a.s., 1^" = VF" + a" and W^"(t) G Us^ (G) for all t > 0, 

(ii) P^-a.s.,W^it)=X-it) + E^^JJ^,,^f'^iW^is-),W^{s))dY,^is)for 

all t>0, where for each i £l, 7*'*^ iM'^ x M'^ — > M"' is Borel measurable and 
U'^{y,x)\\ = lfor allx,yeW', 

(iii) y" = y" + /?", where /3" is locally of bounded variation and P'^-a.s., 
for each i ^Z, 

(a) y«(0)=0, 

(b) Y^ is nondecreasing and AY^{t) < (5" for all t > 0, 

(c) Yp{t) = f^Q^^^ ^{w-^{s)(iUsn{dG,ndG)} 

(iv) 5" — > as n ^ 00, and for each e > 0, there is r]^ > and > 
such that for each i£l, ||7*'"'(?/, x) — 7*(a;)|| < e whenever \\x — y\\ < % and 



18 W. KANG AND R. J. WILLIAMS 

(v) a" — > and V(/3") — > in probability, as oo, 

(vi) {X"} is C-tight. 

Remark. A simple case in which (iv) above holds is where 7*'"(y,x) = 
7*(y). In (v), V(/3") is the total variation process for /?" (cf. Section 1.1). 

The following theorem will play an important role in the proof of the in- 
variance principle. It will be used to show that a sequence of processes sat- 
isfying suitably perturbed versions of the defining conditions for an SRBM 
[cf. (i)-(vi) above] is C-tight. 

Theorem 4.2 (C-tightness). Suppose that Assumption 4.1 holds. Define 
= (Ty"-, X", y") /or each n. Then the sequence of processes {Z^}^^i is 
C-tight. 

Remark. Note that C-tightness of {W}, {X"} and {Y"} implies C- 
tightness of {2"} (see Chapter VI, Corollary 3.33 of [12] for details). 

Proof of Theorem 4.2. References here to (i)-(vi) are to the condi- 
tions in Assumption 4.1. 

Simple algebraic manipulations yield P^-a.s., 



(33) Wit) = X"(t) +y f f '"'(Wis-), Wis)) dYl'is) 

(34) =x«(t)+y"(t)+E / l\W^'{s))dYl^{s), 



where 



(35) X"(t)=X"(i)+f-a"(i) + E / Y'''{W^{s-),W^{s))df5ns)] 



and 

^"W =E / h''''{W^{s-),W^{s))-Y{W-{s)))dY:^is) 

,;gj^{0,i] 

+ E / {Y{W^{s))-j\W^{s)))dYr{s). 



(36) 



The hypotheses on a", the total variation process V(/3"') of Z^", and the 
fact that ||7*''^(y, x) |[ = 1 for all x,y G M"^ and each i € T, imply that the 
process 

-q"(-) + E / y'"(VF"(s-),W^"(s))(i/3f(s) 
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converges to in probability as n — > oo. Combining this with the fact that 
is C-tight, we obtain that is C7-tight. 

Recall the positive nondecreasing function n(-) from Theorem 4.1, and 
the constants a, L and functions R{-) and D{-) from assumptions (A1)-(A5) 
in Section 3. Recall also that Po = 

Fix p,e,r],T such that < p < min{^, ^^^^}, e > 0, > and T > 0. 
By assumption (A3), there is a constant ri > such that 

(37) L»(r) <min||,e| for all r € (0, ri]. 

Since Il{u) — s- as u ^ 0, there are constants < < r2 < min{ri, f , ^} 
such that 

(38) n(r) < y for all r € (0, rg] . 

By (iv), there are < e < min{|, |, -^j and no > such that for all n > no, 

(39) sup sup rnax ||7*'"(?/, x) — 7*(x) II < 



\\y—x\\<2e 



iei 6Ir2 



By (iv)-(vi), and Proposition 1.1, there exist an integer ni > no, a con- 
stant M^^T > and A G (0,T), such that for all n > ni, 

(40) P"( sup ||X"(s)||<M^,t|>1-W2, 

(41) P"{7x;T(X",A)>e}<r?/4, 

(42) P"( sup ||a"(s)||<-^A^j>l-W4, 

lo<s<T 6ILr2 8 J 

(43) J"<min(^,^, / 

^ ^ I 3 21' 8(1 + 1) 21 J 

To prove C-tightness of {W""} and (and hence of {VF"}, {V}), 

by Proposition 1.1, it suffices to show that there exists a constant N^^t > 
such that for all n > ni , 

(44) -P"{'Wt(W^",A) >e} <??, 

(45) P''{wT{Y^,X)>e}<rj, 

(46) P"| sup ||W^"(s)||<iV^,Tj>l-r/, 

Lo<s<T J 

(47) sup ||y"(s)|| <iV^T| > 1 
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For each n > 1, let F" be a set in such that P'^iF'^) = 1 and on F", 
properties (iii)(a)-(c) hold, (33)-(36) hold, and W'^{t) € \Jf,^{G) for all t > 0. 
Fix a t such that < t < t + A < T. Let 

(48) r" = inf {s > t : W'^is) € Usn {dd n dG) for some i G 1]. 
For each n>ni, let 

H^ = {wTiX^,X)<i, sup ||a"(,)||<-!^A^, 
I o<s<T 6ILr2 8 

(49) 

sup ||X"(s) II <M^,T I nF". 

0<s<T ' } 

Then by (40)-(42) and the definition of F", 

(50) P{i?"} > 1 - r/. 

Fix cj" G if". By the definition of wt{x^ A) in (3), we have that, 

(51) sup ^ ||X"(s,cj")-X"(r,cj")|| <e. 

Now there are two cases to consider for n>ni and u, v fixed such that 
t<u<v <t + \. 

Case 1. G {r" > ti}. In this case, by (iii)(c), y"(-,u;") does not increase 
on the interval {u,v], that is, Yl'{v,uj^) — Y^^{u,u;'^) = for ah i G 1. Then 
by (34) and (36), 

(52) W^"(z;,a;") - W^"(u,u;") = X"(z;,u;") - X"(n,a;"). 
Hence, by (51), 

||M^"(?;,w") -Ty"(u,w'^)|| < sup ^ ||X"(s,cj")-X"(r,cj")|| <e<e/8, 

r,sG[i,t+A] 

and we also have 

||y"(7;,cj") -y"(7x,a;")|| = 0<e/2. 

Case 2. G {r" < t;}. Then there is an i G X such that W^{t"- ,lo"-) G 
Us"{dGi ndG), since the set Usr^{dGi HdG) is closed and W^{-,uj"-) is right 
continuous. It follows that there is some xq G dG (which depends on w") 
such that Ty"(r",u;") is in the closed ball Bsn(xo) C Bp{xo). To apply the 
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oscillation inequality in Theorem 4.1, we first prove the following: 

(53) W^"(r, cj") € Bp{xo) for ah r satisfying r" < r < v. 
For the proof of (53), let 

(54) C = mf{r > r" : W^"(r, w") ^ 5p(xo)} A 

By the definition of VF"(r, w") € Bp{xQ) for each r G [r",^"). In order to 
apply the oscillation inequality in Theorem 4.1 on the time interval [t",^"), 
we show that 

(55) Z)(n(Osc(X"(.,a;") + y"(-,^"), [r" ,n) + < f • 

For each r € (0,r], by (i)-(iii) and (33), (49), (43), we have that 
||VF"(r-,u;'^)-VF"(r,cj")|| 

< ||H^"(r-,u;") - W^"(r,a;")|| + ||a"(r-, w") - a"(r, w")|| 

< ||AX"(r,w")||+2 sup ||a"(s)|| +15" 

0<s<T 

^ e e 
<e + - + -<2e. 

Hence by (39), for each r € (0,T], 

(56) ||y."(T^"(r-,u;"),VF"(r,u;"))-y(VF"(r,^"))||<-^. 

6Ir2 

By (36), (56), Assumption (A4), (i) and (49), we have that for any si, S2 
such that u < si < S2 < V, 

\\V^{s2,u;^)-V^{si,u;n\\ 

<IZ/ ||y'''(l^"(r-,u;'^),VF"(r,a;'^)) 

^'"'"'^ -f (l^"(r,w"))||(iy,'^(r,a;") 



+ E / \h\W-{r,u;-))-f{W-{r,u;-mdY:^{r,Lon 



(si,S2] 



(57) 



^E7^(^"(^2,^")-ll"(5i,^")) 



+ E / L||T^"(r,a;")-l^"(r,u;")||dF,"(r,^") 

jgI-'{^l,S2] 

<E7^^(>'"(^2,^")-?r(^i,^")) 
+ E^7^1^(^"(«2,u;")-i;'^(^i,u;")) 



(59) 
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<^\\Y-[s2,u:^^)-Y'\s^,u:^)\\. 

Let 

(58) a" = mf{s>r":Osc(y"(-,u;"),[r",s)) >r2}. 

Note that Osc(y"(-,ti;"), [r", s)) as a function of s defined on (r'^,oo) is 
left continuous with finite right Umits and is nondecr easing. By the right 
continuity of Y^, we know that 

Osc(y'^(-,w"),[r",s))^0 as sir". 

Thus, a" > r", Osc(y"(-, cj"), [r",cj")) < and on {cj" < oo}, Osc(y"(-,u;"), 
[r",(j"])^> ra. By (57), (51), (43), the choice of e, and since t < r" < ^" < 
f < t + A, we have 

Osc(X"(-,u;") + y"(-,cj"),[r",rAcj")) + 5" 

<Osc(X'^(-,cj"),[r",e" Acj'^)) 

+ Osc(y"(-,u;"), [r", ^ A a")) + (5" 

<Osc(X"(-,w"),[r",e" AfT")) 

+ — Osc(y"(-, tj"), [r", r A (t")) + (5" 

<e + :^r2 + <5"<r3. 
3r2 

Then by (38) and the monotonicity of -D(-), we have 

Z)(n(Osc(X"(-, cj") + y"(-, cj"), [r", C A cj")) + 5")) 

(60) 

<i^(^^j<Z)(r2)<Z?(ri)<|. 

We claim that 

(61) fx" > . 

To prove (61), we proceed by contradiction and suppose that a" < Then 
by (60), with x = X"(-,u;") + y"(-,a;") and <5 = (5", condition (iv) of The- 
orem 4.1 holds with [s,t] = [r"',c7"' — l/m] for all m sufficiently large. By 
applying Theorem 4.1 and letting m^oo, we obtain using (34), (38) and 

(59) that, 

Osc(y"(.,u;"),[r",a")) 

(62) < n(Osc(X"(-, a;") + F"(-, w"), [r", C A d")) + <5") 

<n(r3)<^. 
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By (62), (iii)(b) and (43), we obtain that 

Osc(r«(-,^"), [r",a"]) < ^ + 15" < T2. 

This contradicts the fact that Osc(y"(-,u;"), [r", cr"]) > r2 on {d" < oo}, 
and so (61) holds and (55) fohows by (60). 

By applying Theorem 4.1 on [r'",^" — 1/m] and then letting m — > oo, we 
obtain using (61), (59) and (38), that 

Osc(W'^(-,w"),[r",e")) 

< n(Osc(X"(-, a;'^) + F'^(-, w"), [r'^, A a'')) + S'') 

2 ' 

and similarly, 

(63) Osc(y"(-,u;"),[r",n)<y- 
Then we have 

\\W'-{£r-,uj^)-xo\\ 

< WW'^ie-,^! - W^"(r",u;")|| + ||W^"(r",cu") -xoll 

<!^ + 5r 

Using hypotheses (ii), (iii)(b), and (33), (51), we obtain 

iiH?'^(r,w") -w^"(r-,'^'')ii 

< ||X"(e",cj")-X"(C-,cj")|| 

+5]iiy'"(p^"(r-,^"),vF"(r,u;"))ii 

<e + I(5". 



Hence 



|VF'^(e'^,w")-xo|| < \\W''{C-,^'^)-xo\\ 

+ \\W''{C,uj'')-W''{C-,(^'' 

<,-+(I + l)5" + | 
<p/8 + p/8 + p/8<p/2. 
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It follows from this that ^" = f and (53) holds, as desired. 
Then, by (33), (51), (iii)(b), (iii)(c), (63) and (43), we have 

||H?'^(t;,u;")-Ty"(n,w")|| 

< sup |lX"(.,u;")-X"(r,u;")||+^(y,"(^,u;")-y,"(u,a;")) 

r,s&[u,v] jgj 

<e + ^(yi"(?;,w") - yi"(nVr",u;")) 
iei 

(64) +^(y,"(nVr'^,u;")-y,"(t/,a;")) 

< e + I Osc(y"(-, u;"), [u V r", v)) + ^ Ayi"(i;, w") + 

<e + I^ + I(5" + I,5"<- + - + — + — = - 
- 2 8 4 16 16 2 

and 

||y"(t;,u;")-y"(n,a;")||<5](y,"(7;,.;")-y,"(n,a;")) 

<^(yr(r;,a.")-yr(nVr",u;")) 

+ ^(y,"(n V t",u;") - yi"(n,a;")) 
e 

<2- 

Here we have used the fact that y does not increase on (u, r" V u) and can 
jump at most by 6"" at r", by the definition of r" and (iii)(c). 

On combining the results from Case 1 and Case 2, we obtain that for each 
n > ni, 

(66) supj sup ^ |lW^"(t;,u;")-Ty"(u,a;")||:0<i<t + A<rj<^<e 



(65) 



and 

(67) supj sup ^ ||y"(7;,cu'")-y"(u,a;")|l:0<t<t + A<r|<|<e. 

Hence since G if" was arbitrary, by (50), we have that (44) and (45) hold 
for all n > ?ii . 

Next we show that there is a constant N^j^t > such that (46) and (47) 
hold for all n>ni. By (66)-(67) above, we have that for each n>ni, oo^ ^ 
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F", t such that < t < t + A < T and t<u<v <t + X, 

(68) ||M^"('y,u;")-W^"(t/,a;")|| <e 
and 

(69) ||y"(7;,w")-y"(u,a;")|| <e. 

Then, for each < s < T, by (68), (69), (49) and (33), we have 
\\W''{s,uj'')\\ < ||M?"(s,w") - i^"(0,cu")|| + ||W^"(0,u;")|| 

< ^ ||t^"(iAAs,w")-T^"((i-l)AAs,w")|[ + ||X"(0,w")|| 

<([T/A] + l)e + M^,T 

and 

||y'^(s,u;")|| < ||f ''(s,u;") - y"(0,u;")|| 
[r/A]+i 

< ll>'"(^AAs,cj")-y"((«-l)AAs,u;")|| 

i=l 

<([r/A] + i)e. 

Here [T/X] is the greatest integer less than or equal to T/X. Let iV^,T = 
([T/A] + l)e + M^,T- Then we obtain that for n>ni and € i^", 

(70) sup ||W^"(s,cj")|| <iV^,r 

0<s<T 

and 

(71) sup \\Y^{s,L0'')\\<Nr,,T. 

0<s<T 

Then by (50), we have that (46) and (47) hold for all n > ni. 

Finally by applying Proposition 1.1, we have the C-tightness of {Ty"} 
and {y"}^It then follows that {(W^", X", y")}^=i is C-tight. Since 2:" = 
(W^",X",y") + (a",0,/3") where a", V(/3") ^ in probability as n ^ oo, 
then is also C-tight. □ 

4.3. Invariance principle for SRBMs. The main theorem of the paper is 
the following. 

Theorem 4.3 (Invariance principle for SRBMs). Suppose that Assump- 
tion 4.1 holds. Define = {W'^,X'^,Y^) for each n. Then the sequence of 
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processes {Z'^}^^i is C-tight and any (weak) limit point of this sequence is 
of the form Z = {W^X,Y) where continuous d- dimensional processes W,X 
and a continuous I-dimensional process Y are defined on some probability 
space {^i,J-,P) such that conditions (i), (ii) and (iv) of Definition 2.1 hold 
with J^t = ■.0<s<t}, t>0. 

If, in addition, the following conditions (vi)' and (vii) hold, then any 
weak limit point of the sequence {Z^''}^^i is an extended SRBM associated 
with the data (G, /x, T, {7*, z G X}, z/). If furthermore the following condition 
(viii) holds, then as n — s- cxd where W is an SRBM associated with 

{G,fi,T,{f,iGl},u). 

(vi) ' {^"} converges in distribution to a d-dimensional Brownian mo- 
tion with drift /i, covariance matrix T and initial distribution u. 

(vii) For each (weak) limit point Z = {W,X,Y) of {Z^'}"^^^, {X{t) - 
X(0) — fit, Tt, t > 0} is a martingale. 

(viii) If a process W satisfies the properties in Definition 2.1, the law 
of W is unique, that is, the law of an SRBM associated with the data 
{G, fi, r, {Y,i E I}, v) is unique. 

Remark. We note that (vi)' implies that (vi) of Assumption 4.1 holds. 

Proof of Theorem 4.3. By Theorem 4.2, we have that the sequence 
{Z''}n=i is C-tight. Let Z = {W,X,Y) be a (weak) limit point of {^"}^=i, 
that is, there is a subsequence {n^} of {n} such that Z"^*" =^ 2^ as A; — 00. It 
also follows that Z"" = (T^"fe,X"sy"'=) ^ Z as /c ^ 00. By the C-tightness 
of {Z^}, we obtain that Z has continuous paths a.s. For the purpose of 
verifying that Z satisfies the listed properties in Definition 2.1, one may 
invoke the Skorokhod representation theorem to assume, without loss of 
generality, that Z'^'' and Z"* converge u.o.c. to Z A; ^00 and V (/?"*) 

converges u.o.c. to k ^ 00. With this simplification, it is easily 

verified that the properties of j^"*^} and {2^"*=} imply that Z has properties 
(ii) and (iv)(a)-(b) of Definition 2.1. For the verification of property (i) of 
Definition 2.1, note that for each k, a.s. for each t > 0, 

VF"fc(t) =X"'=(i) + y" / y'"*(Ty"Hs-),W^"*(s))'i/3r(«) 

+ Y1 I {l"''"'{W^{s-),W''^{s))--f\W''^{s)))dY^'{s) 

The sum of the first two terms on the right-hand side of the above equality 
converges a.s. to X{t) as k ^ 00. The third term on the right-hand side 
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converges a.s. to as A; ^ oo, by property (iv) and the fact that a.s., 
sup - 

sG(0,i] 

< sup ||AX"'^(s)|| +1 sup ||Ay"'=(s)|| ^0 asfc^oo. 

se(0,t] sG{0,t] 

It remains to show that for each f G X and t>0, a.s., 

f f(w''^is))dY^''''{s)^[ Y{W{s))dYi{s) as fc ^ oo. 

J(0,t] J{0,t] 
This fohows directly from Lemma A. 4. 

For the verification of property (iv)(c) of Definition 2.1, it suffices to show 
that for each z E T, m = 1, 2, . . . , a.s. for each t>0, 

(72) Yi{t)= f fmiW{s))dY,{s), 

J{0,t] 

where {fm}m=i is a sequence of real valued continuous functions defined 
on Mf^ such that for each m, the range of fm is [0,1], /m(a^) = 1 for x G 
Ui/raidGi n dG) and /„^(x) = for x ^ U2/rn{dGi n dG). The existence 
of such a sequence of continuous functions {fm}m=i can be shown using 
Urysohn's lemma (cf. [8], page 122). Then (72) is a consequence of Lemma 
A. 4, property (iii) of Y^"* and the fact that 5"*= —>■ as k ^ oo. Indeed, a.s., 
for each t>0, 

YJt) = limy"''(t)= lim / 1,~„. . , ,ar^ r.ar^\^ dYj^'' (s) 

= lim / U{W^''{s))dYl"={s) 

= / U{W{s))dYi{s). 
J(o,t] 

Thus, Z satisfies properties (i), (ii) and (iv) of Definition 2.1 with J^t = 
a{Z{s):0<s<t}, t>0. 

Assuming properties (vi)' and (vii) hold, Z satisfies (iii) of Definition 
2.1. Then Z is an extended SRBM associated with the data (G, fj,, F, £ 
I}, v). If in addition, property (viii) holds, then the law of W is unique. Since 
each weak limit W is an SRBM associated with the data (G,^,F,{7*,i G 
X}, v) and the law of such an SRBM is unique, then by a standard argument, 
VF" =^ VF as n ^ oo where W is an SRBM associated with (G,/i,F, {7'',i G 
□ 

Some sufficient conditions for (vii) to hold are given in Proposition 4.2 of 
[15] for a simpler setting where G = W^^. Two of those conditions generalize 
to our setting here and can be proved in the same manner as in [15]. For 
completeness, we state the ensuing result here. 
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Proposition 4.1. Suppose that Assumption 4.1 and (vi)' of Theorem 
4.3 hold. If, in addition, one of the following conditions (I)-(II) holds, then 
condition (vii) of Theorem 4.3 is satisfied, and any weak limit point of 
{Z"'}'^^i is an extended SRBM associated with (G, /i, F, {7*, i gX},i/). 

(I) For any triple of d-dimensional {Tt}-o-dapted processes {W,X,Y) 
defined on some filtered probability space ($7, .F, {^(}, P) and satisfying con- 
ditions (i), (ii) and (iv) of Definition 2.1 together with the condition that 
X , under P, is a d-dimensional Brownian motion with drift vector /i, co- 
variance matrix T and initial distribution v, the pair {W,Y) is adapted to 
the filtration generated by X and the P-null sets. 

(II) = X" + e^, y" = + e^, W'' = W'' + e^, where e1,e^,el^ are 
processes converging to in probability as n — > 00, and: 

(a) {X^{t) — X^{0)}'i^^i ^-^ uniformly integrable for each t>0, 

(b) there is a sequence of constants in such that 
lim.„^oo = Ai; 

(c) for each n, {X'^{t) — X"'(0) — n'^t,t >0} is a P"- -martingale with 
respect to the filtration generated by (W"^ , X'^ ,Y") . 

In the rest of this work, we focus on apphcations of the invariance prin- 
ciple and in particular on giving sufficient conditions for property (viii) of 
Theorem 4.3 to hold. 

5. Applications of the invariance principle. In Section 5.1, we prove weak 
existence of SRBMs associated with data (G, /i, F, {7*, i satisfying 
(Al)-(A5) of Section 3. This is accomphshed by constructing a sequence of 
approximations whose weak limit points are SRBMs. The invariance prin- 
ciple is used to prove the C-tightness of the approximations and that any 
weak limit point is an SRBM. In Sections 5.2 and 5.3, using known results 
on uniqueness in law for SRBMs, we illustrate the invariance principle for 
certain domains and directions of reflection. 

5.1. Weak existence of SRBMs. 

Theorem 5.1. Suppose that assumptions (Al)-(A5) of Section 3 hold. 
Then there exists an SRBM associated with the data (G, /i, F, {7*, z gZ},^). 

Proof. We construct a sequence of approximations to an SRBM and 
use the invariance principle to establish weak convergence along a subse- 
quence to an SRBM. 

In the following we will use R{-) from assumption (A2), L > from as- 
sumption (A4), a > from assumption (A5), and po = Fix e > and 
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< p < min{^, ^("^/^^ }. By assumption (A3), there is a constant ri > 
such that 

Z)(r) < min|^,ej> for all r € (0,ri]. 

Recall the properties of n(-) from Theorem 4.1. Since n(u) — > as ti — s- 0, 
there are constants < < r2 < minjri, f , ^} such that 

n('^) < y for all r G {0,rs]. 

Fix e and 5 such that < e < min{|, |, ^, 24iir2 i ^^"^ < < min{^, ^, 
p ±\ 

8(1+1) ' 21/- 

We will construct a d-dimensional stochastic process and an I-dimensional 
"pushing" process , such that approximately satisfies the conditions 
defining an SRBM for the data (G, F, {7*, i Gl},^) (cf. Assumption 4.1). 
The idea for this construction is to use a Brownian motion X with drift 
vector fi, covariance matrix F and initial distribution v. Away from dG, the 
increments of are determined by those of X. For any time t > such 
that W^{t—) G dG, we add an instantaneous jump to W^{t—) to obtain 
W^{t) G G. Here W^O-) = X{0). The size of the jump is such that W^{t) 
is a strictly positive distance (depending on 6) from the boundary of G. 
The jump vector is obtained as a measurable function of W^{t—). To ensure 
the measur ability, each point x on dG is associated with a nearby point x, 
chosen in a measurable way from a fixed countable set of points in dG. The 
jump vector for x is one associated with x. We now specify the mapping 
X — > 5; and the associated jump vector more precisely. 

By assumption (A5)(ii), for each x G dG, there is c(x) G such that 

(73) ^ Ci{x) = l and min ( ^ Ci{x)j^{x),n^ (x) ) > a. 

By (73), Lemma 2.1 and the fact that n*(-) is continuous on dGi for each 

1 £l, we have that for each x G dG there is G (0, 5) such that for each 
yeBr,{x)ndG, 

(74) i{y)cl{x) 
and 

(75) niin ( ^ Ci{x)f{x),n^{y))>^. 

It follows, using the G^ nature of dGi and the fact that n*(y) is the inward 
unit normal to dGi at y G dG for each i G I{y), that (by choosing r^; even 
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smaller if necessary) for each x G dG there is m{x) > and Vx € (0, 6) such 
that for each y G Br^x) n dG, (74)-(75) hold and 

(76) y + \ Ci{x)-f\x) G G for ah A G (0, m{x)). 

Let B°^{x) denote the interior of the closed ball Br^{x) for each x G dG. The 
collection {-B°^ (x) : x G 9G} is an open cover of dG and it follows that there 
is a countable set {xk} such that dG C Ufc-Srj;^. i^k) and {x^} n -BAr(O) is a 
finite set for each integer A'^ > 1. We can further choose the set {x^} to be 
minimal in the sense that for each strict subset C of {xk}, {Br^{x) :x G G} 
does not cover dG. Let = {Br^^{xk) \ i{Ji=i Br^-ixi)) fl dG for each k. 
Then ^ for each k, {Dj.} is a partition of dG, and for each a; G dG 
there is a unique index i{x) such that x G -Dj^^;)- For each x G M'^, let 

c, if X ^ dG, 

Ci{x), iixedG. 

Note that for all x G M*^, 

(77) \\x-x\\<6. 
For each i G X and x G M'^, let 

(78) y'^(x)=f(x). 

The mapping x — > x is Borel measurable on M'^ and hence Y'^ is a Borel 
measurable function from into M"^. 

We construct {W^,Y^) as follows. Let X defined on some filtered proba- 
bility space {i},J-',{J-t},P) be a d-dimensional {^j}-Brownian motion with 
drift n and covariance matrix F such that X is continuous surely and X{0) 
has distribution v. Let 

Ti = mi{t>0:X{t)edG} 

and 

W^{t)=X{t), Y^{t) = forO<t<ri. 

Note that W^{ti—) exists on {ri < oo} since X has continuous paths and 
in the case that n = 0, W\0-)=X{0). On {n < oo}, define 



0, i^liW^in-)), 



and 

H'-'(ti) = X(ti) 



SeJ(VF'5(Ti-)) 
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So have been defined on [0,ri) and at ti on {ti < oo}, such that: 

(i) W'{t) = X{t) + E.gxf •'(W'(0-))y/(0) + E.gx/(o,i] 
Y'\w\s-))dYl{s) for all t G [0,ri] n [0,oo), where W\0-) = X{0), 

(ii) W^t) eG ior t(£[0,n]n[0,oo), 

(iii) for i £l, 

(a) y/(o)>o, 

(b) y/ is nondecreasing on [0,ri] n [0, oo), 

(c) Yfit) = y/(0) + /(o^,] l{H^^(.)ef/,,(aG,naG)} for t G [0, n] n 

[0,oo), 

(iv) \\AY\t)\\ = \\Y^{t)-Y^{t-)\\ <5fortG [0,ri]n[0, oo), where y^(0-) = 

0. 

Note that (iii)(c) above contains the expression W^{s) G U2s{dGindG). The 
reader may wonder why 25 appears instead of 6. The reason is that at a jump 
time s of y/, W^{s-) G dGi n dG and so 

dmi{W\s),dGir\dG) < \\W^[s) - W^[s-)\\ + \\W\s-) -W\s-)\\ 

< 5 + 5 = 26. 

Proceeding by induction, we assume that for some n > 2, ri < • ■ • < t„_i 
have been defined, and ,Y^ have been defined on [0,r„„i) and at r„_i 
on {r„_i < oo}, such that (i)-(iv) above hold with t„„i in place of ri. Then 
we define r„ = oo on = oo}, and on {Tn-i < oo} we define 

Tn = inf{t > Tn-1 : W\Tn-l) + X{t) - X(t„_i) G OG}. 
For Tn-i <t<Tn, let 

y^(t) = y^(T„_i), 

T^'5(t) = 1^-5 (t„_i) + X{t) - X(r„„i), 
and on {t„ < oo}, let 

fv;:'(T„-i), iinw'(T,-j), 

and 

H/'(t„) = W'{t„-) 
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In this way, , have been defined on [0, r„) and at t„, on {r„ < 00} such 
that (i)-(iv) hold with r„, in place of ri . 

By construction {Tn}^=i is a nondecreasing sequence of stopping times. 
Let r = lim„_>ooTn. On {r = 00}, the construction of {W^ , Y^) is complete. 
We now show that {r < 00} = 0. In fact, if {r < 00} 7^ 0, let u; € {r < 
00}. The above construction gives {W^{-,uj), Y^{-,uj)) on the time interval 
[0,T{io)). For each t S [0,t{uj)), we have 

W'{t,u;)=X{t,u;)+J2j''\w\0-,u;))Yf{0,u;) 
i&i 

(79) 

+ E/ y^'{W'{s-,u;))dYfis,u;). 



Since X is continuous on [0, 00), ||7*'''(a;)|| = 1 for each x E M"^ and ^i^jY^ {{), 
uj) < 5, there are constants A E (0,r(u;)) and M > (depending on uj) such 
that 



(80) Wri^)iX{.,u;)+Y'\W'{0-,u;))YfiO,u;),X)<e 
and 



(81) sup 

0<t<T{uj) 



X{.,u;)+y^Y'\w\0-,u;))YfiO,u;] 



where w.{-, ■) is defined in (3). By the choice of e, 6 made at the beginning of 
this proof, (77)-(78) and the uniform Lipschitz property of the 7*(-)) ^ £ ^T, it 
follows that (39) and (43) hold with Y'^{y) and 25 in place of Y'"'{y,x) and 
5" , respectively. Then by similar pathwise analysis to that used in Case 1 and 
2 of the proof of Theorem 4.2, with W"' = W"- = W\ = 0, 7*''^(y,3;) = 
y,<5(y) for j^ach i G I and x,y(^W^, X"^ = X + Y.i&il''\W\<d-))Yl{{)), 
yn ^ Y'' = Y^ - Y\0), /3" = Y^{0) and -5" = 26, we obtain that (71) 
holds for any T < t{uj) with = uj, Nrj^r = i[T{u;)/X] + l)e + M. It follows 
that supjgjSup^g[o,r(a;)) ^/(S)'^) IS finite. By the nondecreasing property of 
Yf{-,uj) on [0,r(a;)) for each i €T, y/(r(a;) — ,u;) exists and is finite for each 
i G X. Then by (79) and the continuity of X, we see that W^{t{oj) — ,uj) exists 
and is finite. By the construction of Y^ and the fact that J2iex{x) = 1 
for all X G dG, we have that 



00 



(82) ^y/w.)-,.)=i:=2:^^^f^Ai. 

igX n=l 

Since Tn{uj) ] t{uj) as n ^ 00 and {t{uj)—,uj) exists, it follows that 
{W^{Tn{uj) — ,uj)}'^^i converges to W^{t{lo) — ,uj) G dG as n ^ 00. Conse- 
quently, {VF'5(r„(cj)-,cj)}^= ]^ is a bounded sequence in dG and so by the 
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definition of the sets {Dk} which form a partition of dG, there is a finite set 
C such that 

{Ty^(T„(a;)-,a;)}~iC (J Dj,. 

k£C 

Hence, 

(83) inf m{W^{Tn{oj)-,uj)) < inf m{xk) > 0, 

n>l keC 

and so the right-hand side of (82) is infinite. On the other hand, since 
suPjgjSup^g[o,T(aj)) ^/(•S)'^) is finite, the left-hand side of (82) is finite. This 
yields the desired contradiction and so {r < oo} = and we have con- 
structed {W^Y^) on [0,oo). 

From the construction above, we can see that and are well-defined 
stochastic processes with sample paths in D([0, oo), M'^) and Z)([0, oo),M^). 
They are adapted to the filtration generated by X and satisfy (i)-(iv) above 
with [0,00) in place of [0,ri]. 

Consider a sequence of sufficiently small 5's, denoted by {6"'}, such that 
(5" J. as n — > 00. For each 6"', let {W^" ,Y^") be the pair constructed as 
above for the same process X. By the above properties and the fact that for 
each i £2 and x, y G M'^, 

(y) _ y (a;)|| < (y) _ y (^) II < L\\y - x\\ < L{5^ + \\y- x||), 

we obtain that Assumption 4.1 holds with = = W^" , a" = 0, Y'^iu^ = 
^i,S" (y) for gj^ci^ i G Z and X, y G M"', X'^ = X + J2iei f'^" (W^" i^-Wf" (0), 
yn ^ y5" ^ yn ^ y5" _ y<5" (Q) , = y^" (Q) and 25*^ in place of 5" . By invok- 
ing the first part of Theorem 4.3, we obtain that {Z^"}^^^ = {{W^",X^", 
y^")}^]^ is C-tight and any weak limit point Z of this sequence satisfies 
conditions (i), (ii) and (iv) of Definition 2.1 with J^t = ^"{^(s) - < s <t}, 
t>0. Note that condition (vi)' of Theorem 4.3 holds trivially. Furthermore, 
M^" = {X^"{t) - X^"{0) - /it, t > 0} = {X{t) - X{0) - fit,t> 0} is a mar- 
tingale with respect to the filtration generated by X. Since W^" , Y^" are 
adapted to this filtration, it follows that M^" is a martingale with respect 
to the filtration generated by W^" , X^" , Y^" (which in fact is the same as 
that generated by X). For each t > 0, X^"(t) - X^"(0) = X{t) - X{0) and 
so trivially this forms a uniformly integrable sequence as n varies. It fol- 
lows from Proposition 4.1 that condition (vii) of Theorem 4.3 holds. Hence, 
any weak limit point of {Z^"}^^i is an extended SRBM with the data 

(G,/i,r,{y,ieT},i/). □ 

5.2. SRBMs in convex polyhedrons with constant reflection fields. Exis- 
tence and uniqueness in law for SRBMs living in convex polyhedrons with a 
constant reflection fleld on each boundary face has been studied by Dai and 
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Williams [4]. In this subsection, we state a consequence of our invariance 
principle using the results in [4] to establish uniqueness in law. In this case, 
G is defined in terms of I (I > 1) d-dimensional unit vectors {n^,i € 1} and 
an I-dimensional vector (3 = {f3i , . . . , such that 

(84) G = {xeR'^: {n\ x) > /3, for all i G I}. 

It is assumed that G is nonempty and that the set {(n^,/3i), . . . , (n^,/3i)} is 
minimal in the sense that no proper subset defines G. For each i £l, let Fi 
denote the boundary face: {j; G G: = Pi}. Then, n* is the inward unit 

normal to -Fj. A constant vector field 7* of unit length specifies the direction 
of reflection associated with Fi . 

Definition 5.2. For each 0^}CcI, define Fic = Clinic ^i- Let F0 = G. 
A set /C C T is maximal if /C 7^ , Fic 7^ and F)c ^ Fj^ for any AT D /C such 
that KL^K. 

In [4], Dai and Williams introduced the following assumption. 
Assumption 5.1. For each maximal KLcl, 

(S.a) there is a positive linear combination n = X^igx: (bi > \/i (z IC) 
of the {n^,i G /C} such that {n, 7*) > for all i G /C, 

(S.b) there is a positive linear combination 7 = X^jgC (^Cj > Vi G /C) 
of the {7*,i G /C} such that (n*,7) > for all i £ IC. 

Remark. For the given G and constant vector fields {7^,1 GT}, As- 
sumption 5.1 is equivalent to assumption (A5). 

Definition 5.3. The convex polyhedron G is simple if for each )C Cl 
such that /C 7^ and F/c 7^ 0, exactly |/C| distinct faces contain Fic- 

Remark. The polyhedron G is simple if and only if JC is maximal for 
every /C such that ^ fC Cl and Fic 7^ 0- It is shown in [4] that when G 
is simple, (S.a) holds for all maximal /C C X if and only if (S.b) holds for all 
maximal /C C X. 

Dai and Williams [4] showed that Assumption 5.1 is sufficient for exis- 
tence and uniqueness in law of SRBMs living in G with the reflection fields 
{7*,i G X} and fixed starting point. [They also showed that condition (S.b) 
holding for all maximal /C C X is necessary for existence of an SRBM starting 
from each point in G. Consequently, when G is simple. Assumption 5.1 is 
necessary and sufficient for existence of an SRBM starting from each point 
in G.] This yields the following consequence of our invariance principle. 
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Theorem 5.4. Let G be a nonempty domain such that G is a convex 
polyhedron of the form (84) (with minimal description), and let {7*,^ €X} be 
a family of constant vector fields of unit length satisfying Assumption 5.1. 
Suppose that Assumption 4.1 and (vi)', (vii) of Theorem 4.3 hold. Then 
=^W as n —5- 00 where W is an SRBM associated with (G,/i,r, {7*,? € 

Proof. Clearly (Al) holds. Assumptions (A2)-(A3) hold by Lemma 
A. 3. Since for each i G T, 7*(-) is a constant vector field of unit length, as- 
sumption (A4) holds trivially. Assumption (A5) is implied by Assumption 
5.1. Hence by Theorem 4.3, the only thing that we have to check is condition 
(viii) of Theorem 4.3, that is, uniqueness in law for SRBMs in convex poly- 
hedrons with constant reflection flelds of unit length. But this is proved in 
Theorem 1.3 of [4] for a fixed starting point in G and follows by a standard 
conditioning argument for the initial distribution v. □ 

5.3. SRBMs in bounded domains with piecewise smooth boundaries. Dupuis 
and Ishii [6] have established sufficient conditions for the existence and path- 
wise uniqueness of reflecting diffusions living in the closures of bounded 
domains with piecewise smooth boundaries. In this subsection, we state a 
consequence of our invariance principle using the results in [6] to establish 
uniqueness in law. 

Theorem 5.5. Let G be a bounded domain and {7',? S X} be a family 
of reflection fields that satisfy assumptions (Al)-(A4) and (A5)' in Section 
3. We further assume that for each i S I, 7*(') is once continuously differ- 
entiable with locally Lipschitz continuous first partial derivatives. Suppose 
that Assumption 4.1 and (vi)', (vii) of Theorem 4.3 hold. Then ^ W 
as n —i- 00 where W is an SRBM associated with (G, /i, F, {7*, i € T}, ;/). 

Remark. We remind the reader that in view of Lemma 3.1, to verify 
condition (A5)', one only needs to show that (i) or (ii) holds for all x € dG. 
However, as can be seen from the proof below, both forms of the condition 
can be useful. 

Proof of Theorem 5.5. This theorem follows from Theorem 4.3 and 
uniqueness in law for the associated SRBMs. The latter follows by a standard 
argument from the pathwise uniqueness established in Corollary 5.2 of [6] for 
their Case 2. The conditions required for that case are satisfied in particular 
because (A5)'(ii) implies condition (3.8) of [6]. That condition (3.8) readily 
implies condition (3.6) of [6]; and, by [5], under the additional smoothness 
assumptions imposed on the 7* in the statement of our theorem, condition 
(3.8) also implies condition (3.7) in [6]. In addition, (A5)'(i) implies that 
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for each x E dG, {Y (x) , n'^ (x)) > for each i S and furthermore, since 

(A5)' imphes (A5), we have by (A5)(i) that the origin does not belong to 
the convex huU of the {'y^{x) :i G I{x)}. □ 

APPENDIX: AUXILIARY LEMMAS 

Lemma A.l. Suppose that G is hounded. If assumption (Al) holds, then 
assumption (A2) holds. 

Proof. To see this, suppose G is bounded and assumption (Al) holds. 
Fix e € (0, 1). For each i^X and z G dGi n dG, by the property of dGi, 
there is a neighborhood Vz of z and a constant R{e,i,z) > such that for 
all X G Vz n dGi n dG and y £ Gi such that ||x — y|| < R{e, i, z), 

(85) {n^{x),y — x) > —e\\y — x\\. 

Assumption (A2) then follows by a standard compactness argument. □ 

Lemma A. 2. Suppose that G is a nonempty bounded domain satisfying 
(5), where for each i &I, Gi is a nonempty domain. Then assumption (A3) 
holds. 

Proof. We prove the lemma by contradiction. Suppose that assumption 
(A3) does not hold. Then, since there are only finite many J' Cl, J ^ 0, 
there is an e > 0, a nonempty set J' Cl, a sequence {r„} C (0, oo) with r„ — > 
as n ^ oo, a sequence {xn} C K'^ such that for each n, Xn G Cljej {dGj PI 
dG) and dist{xn, Hj^jidGj (1 dG)) > e. But since G is bounded, {xn} is 
bounded and without loss of generality we may assume that 
n—)-QO for some x G W^. It follows that x G CljejidGj Ci dG), since for each 

dist(x, dGj n dG) < — + disi{xn, dGj n dG) < — x|| + r„ ^ 

as n ^ oo. This is inconsistent with Xn^ x and dist(x„, Cljeji^Gj fl dG)) > 
e. □ 

Lemma A. 3. Suppose (Al) holds where 

(86) Gi = {xGM'^:(n*,x) >/3i} fori el, 

{n^,i G X} is a finite collection of d- dimensional vectors of unit length, and 
for I = \I\, P = (/?!, . . . , /3i)' is an I-dimensional vector. (Thus, G is a convex 
polyhedron.) Assume that for each i el, dGi H dG ^ 0. Then assumptions 
(A2) and (A3) hold. 
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Proof. Assumption (A2) holds automatically since G is convex. In or- 
der to show that assumption (A3) holds, we just need to show that for each 
J^CT with ^7^0, 

(87) supJ dist [ X, Pi {dGj n aC) J : x G fj C/^ {dCj n aC) I ^ 

as r ^ 0. Fix J C T such that J ^0. Then njgj(<9Gj n^G) is the collection 
of all solutions x E M'^ to the following system of linear inequalities: 

(n* ,x)>(3i for all i G T, 

(LS) 

{-n\x)>-(3i for ah i € J. 

Suppose that C\j^j{dGj CidG) / 0, that is, (LS) has at least one solution. 
By a theorem of Hoffman [11], with supporting lemmas proved by Agmon 
[1], there is a constant C > (depending only on {n^,i S 1} and not on /3) 
such that for any x € M'^ there exists a solution xq G M'^ of (LS) with 

(88) ||x - xoll < C [ ^(A - {n\ x))+ + ^ (-ft - {-n\ x))A . 

For r > 0, any x G Pljej Ur{dGj n dG) satisfies the following: 
(n*, x) > (3i — r for all i G X, 

(r-LS) 

{-n\ x) > -Pi - r for all i G J. 
Then by (88), there is xq G Cljf^jidGj n ^G) such that 

distfx, fl {dGjOdGU < \\x - xo\\ <2C\I\r. 

It follows that (87) holds when Cljeji^Gj n dG) 7^ 0. 

Now suppose that HjejidGj CidG) = 0, that is, (LS) has no solution. We 
shall use an argument by contradiction to show that Hjej UridGjCidG) = 
for all r sufficiently small. Suppose that this is not true. Then we have that 
r\jfzjUr{dGj n dG) / for all r G (0, 00). As we have seen before, any 
X G r\j£j Ur{dGj n dG) is a solution to (r-LS). We now construct a Cauchy 
sequence. Let xi G f]ji^jUi/2{dGj DdG). Then xi is a solution to (^-LS). 
Since (^-LS) has at least one solution, by the theorem of Hoffman [11] (using 
the fact that the constant C depends only on {n*,i GT}), we conclude that 
there is a solution X2 to (^-LS) such that ||xi — X2II < where C = 2C\2\. 
Continuing in this manner, we can obtain a sequence {x^lj^i such that for 
each n > 1, ||x„ — Xn+i|| < Xn+i is a solution of (^^-LS). The 
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sequence {xn}'^=i is Cauchy. Hence, there is an x* € M"^ such that rE„ x* 
as n — > CO, and x* is a solution to (LS). This contradicts the supposition 
that njGj(9G'j n dG) = 0. Thus we have that fljej Ur{dGj n dG) = for 
ah r sufficiently small, and for such r. 



supJ dist [ X, fl {dGj n aC) ] : X G f| Ur{dGj n dG) 



0, 



by convention. 

Combining the above we see that for each J <ZT with JT" 7^ 0, (87) holds 
and hence assumption (A3) holds. □ 

Remark. In fact, under the assumptions of Lemma A. 3, there is a con- 
stant C > such that D{u) < Gu for each n > and D{-) defined as in 
assumption (A3). 

Lemma A. 4. Given T > 0, functions 0, J^;^ in D([0, 00), M"^), and 
Xiix"}^! -C([0, c)o),M), suppose that supo<s<T ||</>"'(s) — (/)(s)|| — and 
supo<s<T \x"'{s) — x{s) \ — > as n—f 00. Assume that x" is nondecreasing for 
each n. Then for any sequence of real valued continuous functions 
defined on M*^ such that /" converges uniformly on each compact set to a 
continuous function f : M"^ M, we have 



nr{s))dx''{s) 

l{0,t] 

uniformly for t G [0, T] . 



fi4>is))dxis) as 00, 



Proof. By replacing x^'i) and x(-) by x"(-) - x"(0) and x(-) - x(0), 
respectively, we may assume that x"'(0) = x(0) = 0- It is straightforward to 
see by the uniform convergence of {x"} to x on [0,T] that x inherits the 
nondecreasing property of the {x"}- 

By the triangle inequality. 



sup 

0<t<T 



r(0"(5))dx"(5)- / fWs))dx{ 



(90) 



< sup 

0<t<T 



{nr{s))-f{<i>mdx''{s] 



sup 

0<t<T 



/(0(s))(i(x"(s)-x(s)) 



For the first term on the right-hand side of the above inequality, we have 



sup 

0<t<T 



(r(<^"(«)) -/(</>(«))) dx"(s) 



< sup |r(0"(5))-/(0(s))|x"(T), 

0<s<T 
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where the right-hand side member above tends to zero as n ^ oo by the 
uniform convergence of (j)^ to (/> on [0, T] (which imphes uniform boundedness 
of {(/>"} on [0,r]), the uniform convergence of /" to / on compact sets, the 
continuity of /, and the convergence of x"(^) to x{T)- For the second term, 
note that since f{4>{-)) € L'([0, cxo), M), by Theorem 3.5.6, Proposition 3.5.3 
and Remark 3.5.4 of [7], there is a sequence of step functions {z'^}^^ of the 
form 



(91) 



' ^ L 2 ' Z + 1 / 



1=1 



where 1 < < oo, = < t2 < • • • < *f^+i < c» and supQ<g<^ |/(0(s)) 
z^{s)\ — > as k ^ oo. Then 



sup 

0<t<T 



f(4>(s))dix''{s)-x{s)) 



< sup 

0<i<T 



(/(0(.)) -.'=(«)) d(x"(.)-x(^)) 



+ sup 

0<t<T 



z\s)dix''{s)-x{s)) 



< sup |/(0(s))-z'=(s)|(x"(T) + x(T)) 

0<s<T 



Ik 



+ sup J2\z'^(t^)\\ix--xml,At)-)-{x''-Xm'iAt)-)\. 



0<t<T 



i=l 



For fixed k, the last term above can be made as small as we like for all n 
sufficiently large since x" ~^ X uniformly on [0, T] . The desired result follows. 

□ 



Remark. The proof of Lemma A. 4 is a modification of the proof of the 
related Lemma 2.4 in [4]. The difference in assumptions is that in [4] it is 
assumed that (p^ ^ cj) in the Ji -topology rather than uniformly on [0,T], 
x", X ^ ^"([0, oo),M+) rather than x", X S -D([0, oo), M), and there is a single 
function / rather than a sequence {/"'}■ 
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